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We discuss the detectabihty of gravitational waves with a time dependent mass contribution, by 
, means of the stochastic gravitational wave observations. Such a mass term typically arises in the 

cosmological solutions of massive gravity theories. We conduct the analysis based on a general 
c 2 ^ ■ quadratic action, and thus the results apply universally to any massive gravity theories in which 

04 \ modification of general relativity appears primarily in the tensor modes. The primary manifestation 

of the modification in the gravitational wave spectrum is a sharp peak. The position and height of 
r—{ • the peak carry information on the present value of the mass term, as well as the duration of the 

infiationary stage. We also discuss the detectabihty of such a gravitational wave signal using the 
future-planned gravitational wave observatories. 



■ PACS numbers: 



CLi I. INTRODUCTION 



Since the pioneering model of massive gravity was proposed by Fierz and Pauli [l| , numerous attempts have been 
made to introduce a non-zero mass to graviton. This issue has been attracting a great deal of interest, partly 
because graviton mass may provide an alternative explanation for the acceleration of our universe. Namely, instead 
of attributing it to the existence of dark energy whose origin is still unknown, the acceleration may be simply due to 
] modification of gravity. 

• Until recently, however, it has been thought that such modification by non-zero graviton mass is extremely difficult, 
l/^ I typically leading to the emergence of ghost degrees of freedom and various pathologies. As first shown by Boulware 
and Deser the source of these difficulties is associated with a helicity-0 mode in the gravity sector, which is absent 
■ at the quadratic order of the Fierz-Pauli massive gravity but revives at the higher order and acts as a ghost, dubbed 
fS| . the Boulware- Deser (BD) ghost. A remedy for this difficulty was recently proposed by Refs. 0, in which the 
would-be BD ghost is eliminated by adding higher-order correction terms order by order and by resumming them into 
a three parameter expression. Non-existence of the BP ghost in this theory was discussed and proved in Refs. (5l-[ll|. 
Various solutions in the theory were studied in Refs. [12h26| |. Therefore this nonlinear theory of massive gravity can 



> 

in 



X: 



be considered as a nonlinear completion of the Fierz-Pauli theory, which is the simplest and the oldest linear massive 
gravity theory. Having this in mind, we consider it quite natural to suppose that other massive gravity theories or, 
more generally, Higgs phases of gravity might find their nonlinear com plet ions in the future. 

Such exarnples include, but are not restricted to, ghost condensate |27l . and Lorentz-violating massive gravity 
theories [1^, [s^. In these examples, propagating degrees of freedom do not form a representation of 4-dimensional 
Poincare symmetry even in the exact Minkowski background. Instead, they form representations of 3-dimensional 
rotational and translational symmetries. As a result, physical degrees of freedom are classified into scalar, vector 
and tensor parts, according to their transformation properties under the 3-dimensional spatial rotation. Modes in 
different classes may behave rather differently. For instance, in the case of ghost condensate in Minkowski background, 
propagating degrees of freedom in the gravity sector are two tensor modes and one scalar mode. The two tensor modes 
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behave exactly like those in general relativity (GR) while the scalar mode has an unusual dispersion relation, i.e. non- 
relativistic dispersion relation without a mass gap. On the other hand, in a particular class of models among those 
proposed in [30| , propagating degrees of freedom in the gravity sector are two tensor modes only, and their dispersion 
relation is massive and relativistic. From the point of view of the symmetry, i.e. 3-dimensional spatial rotational 
symmetry, we expect that there should be even wider classes of models in which scalar, vector, and tensor parts 
behave differently. 

In the present paper, for definiteness we mainly focus on a class of massive gravity models in which propagating 
degrees of freedom in gravity sector are only two tensor modes with a massive dispersion relation. In this case, 
the scalar and vector sectors behave exactly like in GR and the only modification emerges in the tensor sector: the 
dispersion relation of gravitational waves acquires an effective mass, which in general can be time-dependent. We 
shall discuss an attempt towards a nonlinear completion of such massive gravity models, but in most part of the 
present paper we shall consider this class of models as a purely phenomenological one to be constrained or probed by 
observations and experiments. We conduct the analysis based on a general quadratic action for tensor modes, and 
thus our results apply not only to the theories with precisely two propagating gravitational modes but also to a more 
general class of theories in which modification of gravity appears primarily in the tensor modes. 

Being different from GR only in the tensor sector, observation of gravitational waves and/or their imprints would 
be the most efficient probe for the model. Therefore, in this paper we address the detectability of the effective mass 
of gravitational waves by means of the stochastic gravitational wave observations. 

The rest of this paper is organized as follows. After describing our model in Sec. [ITl we discuss an attempt towards 
nonlinear completion of the model in Sec. IIIII This consideration at the very least shows that the structure of the 
model considered in the present paper is not forbidden by the symmetry. However, readers who are interested only in 
purely phenomenological aspects can safely skip Sec. IIIII We then derive analytic formulas for the power-spectrum of 
stochastic gravitational waves and discuss how to read off information about the effective mass of gravitational waves 
in Sec. IIVI In Sec.|V]we show some numerical results and make comparison with the analytical results of the previous 
section. Sec. |Vl]is devoted to a summary of the results and discussions. 



II. MODEL DESCRIPTION 



In the present paper we shall consider linear massive gravity models that respect the symmetry of Friedmann- 
Robertson- Walker (FRW) universe, i.e. the spatial homogeneity and isotropy. Metric perturbations are thus decom- 
posed into scalar, vector and tensor parts, according to the transformation properties under the 3-dimensional spatial 
rotation. 

We restrict our considerations to a class of models in which propagating degrees of freedom in the gravity sector 
are two tensor modes only. For tensor-type perturbations on the FRW background, the metric can be written as 



where ^ij is the tensor-type perturbations satisfying 



(1) 



(2) 



Here, fi*^ — (51"^)*-' , = Vf^Dj and Dj is the spatial covariant derivative compatible with the metric fiij. 

By assumption, the quadratic action for the perturbation should be invariant under the spatial rotation and trans- 
lation. Hence we can use the following ingredients to construct the quadratic action describing the dynamics of the 
tensor perturbation; 



Ndt, VTid^x, n,., Di, 



1 . 

TV ' 



(3) 



det and an overdot 



provided that spatial indices are properly contracted to form scalar quantities. Here, f2 
represents derivative with respect to the time coordinate t. Since the time coordinate t is invariant under the spatial 
rotation and translation, we can also use general functions of t and their derivatives. 

We also demand that the action should include only up to second time derivatives of "fij in order to avoid ghosts 
in the linear level and that the system respects the spatial parity invariance. 

Thus, after integration by parts, one can show that the general quadratic action is of the form 



r(2) 



I dtd^xNa''VnF{t) 



1 

7V2 



r-'lij + r- 



(4) 
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where 7*-' = Ti-' 7fci, A = D*Di is the Laplacian operator associated with fly, F(t) and c„(<) (n = 0, 1, - • •) are 
general functions of the time coordinate t. Note that F{t) must be positive definite in order to avoid appearance 
of ghost and strong couphng. We can then set = 1 by a field redefinition, to be more precise, by a conformal 
transformation. Finally, in order to describe low energy phenomena such as cosmological evolution of gravitational 
waves at late time, we truncate the series expansion (the sum over n) at the second order spatial derivatives and 
obtain 



r(2) 



PI 



dtd^xNa^Vn 



■^ri,, + (A - 2K) 7,;, - M^^^(t)7^^7,, 



where we have set = 1 by a field redefinition (or a conformal transformation) and we have defined 

Note that physical meaning of Cg and Mqw arc the sound speed and effective mass of gravitational waves. 
Taking advantages of the background symmetry, it is convenient to expand jij as 



where = ^lijk^k^ , A denotes the hclicity state and Y^^ is the tensor harmonics satisfying 

(A + fc2) ^ 0^ JJ^Y^x ^ ^^JY^^ J ^3^^ Q'in''^Y,l{k,x)Yli' {k',x) = 32TT^Sxy6^{k + k'). 



The equation of motion for 7^^^ is then 



7r 



(5) 



(6) 



(7) 



(8) 



(9) 



where the prime (') denotes derivative with respect to the conformal time 77 defined by drj = Ndt/a. Since the equation 
of motion is identical for both polarizations, we omit the index A hereafter. 



III. AN ATTEMPT TOWARD NONLINEAR COMPLETION 

The nonlinear theory of massive gravity recently proposed by Refs. 0, 3 eliminates the renowned BD ghost by 
construction and thus can be considered as a nonlinear completion of the Fierz-Pauli theory, which is the simplest 
and the oldest among all linear massive gravity theories. Having this in mind, we consider it quite possible that the 
phenomcnological model described in the previous section might also find its nonlinear completions in the future. 

In this section we discuss an attempt towards such a nonlinear completion. In particular, we shall derive the model 
with Cg = 1 precisely, based on a non-trivial background in the nonlinear theory of massive gravity. Unfortunately, 
this construction is purely classical and fails at the quantum level [3l[.^ Nonetheless, at the very least it shows that 
the structure of the model considered in the present paper, with Cg = 1, is not forbidden by symmetry. Readers who 
are interested only in purely phenomenological aspects can safely skip this section. 

The covariant action of the nonlinear massive gravity is constructed out of a 4-dimensional metric and four 
scalar fields (a = 0, 1, 2, 3) called Stiickelberg fields. The Stiickelberg fields enter the action only through the tensor 
/pi/ defined as 

Uu = Iab{^ld^V''d,^\ (10) 

where fabif'^) is a non-degenerate, second-rank symmetric tensor in the field space. The spacetime metric g^i, and 
the tensor /^^ are often called physical metric and fiducial metric, respectively. 



^ Adopting a different approach, Ref. |32|| achieved a similar conclusion. 
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The gravity action is the sum of the Einstein-Hilbert action Ieh,a (including a bare cosmological constant A) and 
the graviton mass term Imass specified below. Adding the matter action Imatter j the total action is 



H~ -^mass\_9 1 f ~^ -^matter 

[Sm^^ctj], (11) 



where 



f.f2 

IehA9^.A = J d^x^{R~ 2A), (12) 
Imass[gt^u,ftj.u] = ^^Ipi'^^l j d!^ x^/^ {C2 + a^C^, + a^^C^) , (13) 
and {cT/} (/ = 1, 2, • • • ) represent generic matter fields. Each contribution in the mass term Imass is constructed as 



C2 = ^([/c]^-[/c^] 



2 ' ''^ 

1 /...i3 



-Cs = ^ [[ICf - 3 [/C] [/C2] + 2 [/C3] 
£4 = ^ ([/C]^ -6 [ICf [IC'] + 3 [/C^] % 8 [/C] [/C3] - 6 [/C^] ) , (14) 
where the square brackets denote trace operation and 

/cj; = <5t'- (vFv)''^ • (15) 

The square-root in this expression is the positive definite matrix defined through 

Note that the fiducial metric is not dynamical but fixed by the theory. In other words, the form of the fiducial metric 
is a part of definition of the theory. Of course, an apparent form of the fiducial metric is subject to changes under 
redefinition and dynamics of the Stiickelberg fields. However, such changes are equivalent to coordinate transforma- 
tions of the metric in the field space and thus do not change geometrical properties such as the curvature tensor in the 
field space. For example, if we demand the maximal symmetry, i.e. Poincare, de Sitter or anti-de Sitter symmetry, in 
the field space then the fiducial metric must be Minkowski, de Sitter or anti-dc Sitter metric, respectively, regardless 
of appearances. 

In the present paper we consider a general fiducial metric of the FRW type, specified by the lapse function n{ip^), 
the scale factor a{ip^) and the spacial curvature constant K: 



= -n2(/)a,^0a./ + a2(/)17„-(/)9,^^9.^^ f],,(/) = 5,, + / '! ^ ^ ■ (17) 

1 - KdlmV'-ip"'' 



Needless to say, this class of fiducial metrics include Minkowski, de Sitter and anti-de Sitter fiducial metrics as special 
cases. For example, a Minkowski fiducial metric can be rewritten as an open FRW form and allows non-trivial 
cosmology 18 1. Note, however, that the Minkowski fiducial metric does not allow flat nor closed FRW cosmology [l7| . 
Similarly, an anti-de Sitter fiducial metric allows non-trivial open FRW cosmology but does not allow flat nor closed 
ones (see Appendix [B|) . On the other hand, non-trivial open, flat, closed FRW cosmologies can be realized with a de 
Sitter fiducial metric (see Appendix IbI and see also Ref. [H, [s^l for related discussions). 
In this section we consider a general FRW background plus perturbations. 



= C + ^''^x'^+t:- + iTT^Sfc^'^ + 0(e3), (18) 

where 



gl^'Jdx^^dx'' = ~N{tfdt^ + a{tf^l,,{x^)dx'dx\ n,,(x^) = S,, + J ^ ^ . (19) 

Here, the perturbation 7r° of the Stiickelberg fields was defined through the exponential map, following [sl]. The 
background equations of motion lead to 



'6K , 1 2H 2K 1 ^ ^ ait) , , 
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where a dot represents derivative with respect to t, H = d/{Na) is the Hubble expansion rate of the background 
physical metric, p and P are the matter energy density and pressure, and 



X± = ■ , A± = -m (1 - A±) 3 - A± + a3(l - A±) . (21) 

as + q;4 " 

(See Ref. [11] or Appendix [B] for derivation.) Note that the background physical metric follows the Fricdmann equation 
with an effective cosmological constant A± and thus the universe exhibits self-acceleration if A± > 0. 

Thanks to the exponential map used to define 7r° in (|18p . the action for cosmological perturbations around the 
self-accelerating FRW background manifestly respects the background symmetry. This allows us to decompose the 
perturbations into scalar, vector and tensor parts, according to transformation properties under spatial diffeomor- 
phism. A priori, because of the inclusion of four Stiickelberg fields, we might expect cosmological perturbations in 
the gravity sector to consist of 6 degrees of freedom: 2 tensor, 2 vector and 2 scalar modes. Actually, since the 
theory is constructed in such a way that the would-be BD ghost is completely excised, one would instead expect that 
cosmological perturbations in the gravity sector should consist of 5 degrees of freedom: 2 tensor, 2 vector and 1 scalar 
modes. Surprisingly enough, the gauge-invariant analysis of the quadratic action in [ssj has shown that the scalar 
and vector degrees have vanishing kinetic terms. Assuming these degrees are non-dynamical, we may integrate them 
out, then find that in the scalar and vector sectors, gauge-invariant variables constructed from metric and matter 
perturbations have exactly the same quadratic action as in general relativity. 

Hence, for linear perturbations around the self-accelerating FRW background, difference from general relativity 
arises only in the tensor sector. Assuming that there is no tensor-type contribution from matter fields aj for simplicity, 
the quadratic action for the tensor sector is given by ([5]) with 

I JL^ 

Cs = l, Mi;n. = ±ir~l)mlX±^l + a3 + al-ai, r = (22) 

Here, H = a/{Na) and Hf = a /{no) are Hubble expansion rates of the background physical metric and the fiducial 
metric, respectively. Note that r and thus Mqw are time-dependent functions determined by the ratio H/ Hf. 

We have derived the model described in theprevious section with Cg = 1. As already stated, the derivation is 
purely classical and fails at the quantum level [31| . For this reason we do not consider this construction as a realistic 
one. Nonetheless, at the very least this shows that the structure of the model is not forbidden by symmetry. It may 
be worthwhile focusing on new models such as Ref. (35| and/or investigating cosmological perturbations around new 
classes of cosmological backgrounds such as those found in Ref. [l9l. [26l. IsgI] . towards a successful nonlinear completion. 

IV. STOCHASTIC GRAVITATIONAL WAVE SPECTRUM 

To clarify the influence of Mqw on the observable quantities, we seek approximate solutions to Eq. ([9]) for a generic 
Maw (t) , by considering super- horizon modes to be frozen and sub- horizon modes to follow the WKB-type evolution. 
We then derive approximate analytic formulas for the energy and power spectra. For simplicity we set 

Cg = 1, (23) 

throughout this and the forthcoming sections except for Sec. IIV E 41 where we give a brief comment on the case of a 
general Cg. 

A. Horizon re-entry 

From the action ([5|) and the equation of motion © with Cg = 1, one can easily read off the dispersion relation as 

a.2 = ^ + A4^, (24) 

where a; is the effective frequency with respect to the proper time r defined by dr — Ndt ~ ad-q. Here, since the 
comoving momentum k appears only through the combination ^/k"^ + 2K, we have renamed this combination as k so 
that K does not appear explicitly. We shall still call k (after this redefinition) comoving momentum. We assume that 
the mass squared M^j^ dominates the right hand side of ([M]) at late time and that it does not lead to instability, i.e. 

dt{a'A'4y^) > 0, M^^ > 0. (25) 
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(See Sec. IVland Appendix IB] for a simple example in which these conditions actually hold.) 

The frequency lo given by the dispersion relation (j24p is relevant and the tensor mode oscillates with this frequency 
if and only if it is sufficiently high compared with the cosmic expansion, i.e. ^ . On the other hand, when 
^ , the equation of motion Q allows a growing solution 7^ oc a and thus 7fc ~ const. At a given moment of 
time, a mode with uP 3> is called sub-horizon and a mode with cj^ <C is called super-horizon, as usual. For a 
given mode with comoving momentum k, the time when ~ is called horizon crossing, again, as usual. 

We suppose that scale-invariant super-horizon gravitational waves are generated in the early universe, e.g. through 
inflation. In the present paper we concentrate on the subsequent evolution of gravitational waves, assuming for 
simplicity that the matter energy density and pressure in the background equation of motion ((20)) satisfy the strong 
energy condition {p-\-3P > 0). With this assumption, a^p decreases as the universe expands. Since the curvature term 
{3K/a'^) is small (less than a few percent of 3H^) at present time, this implies that the curvature term is negligible 
throughout the evolution. Hence, before the effective cosmological constant A± starts dominating the universe,^ 
a^H^ decreases as the universe expands and, under the assumption (j25p . super- horizon modes generated in the early 
universe eventually cross the horizon and become sub-horizon. This type of horizon crossing is called horizon re-entry, 
as usual. 

For a mode with comoving momentum k, the time of horizon re-entry is denoted as t = tk and defined as the 
solution to 

Lu^itk) = H\tk). (26) 

We denote the values of a and H at t = tk as Ok and Hk, respectively. For comparison, we define the time of horizon 
re-entry in the massless case (GR) t'j^^ by 

^ = if2 at t = iGi?^ (27) 
and denote the values of a and H at t ^ t'^^ as and H^^, respectively. 

B. Classification of modes 

As defined above, the modes re-enter the horizon when the effective frequency w given by ([24]) exceeds the Hubble 
expansion rate H . Thus, the discussion of the evolution of the modes can be reduced to the determination of the time 
when either the mass contribution a Mqw or the comoving momentum k exceeds the comoving horizon scale a H. 
Hence, in this subsection, we discuss the evolution of a Mow and a H relative to k. 

We define two special scales which characterize classification of the different modes. First, we introduce the critical 
momentum kc for which, the horizon crossing occurs at time t = tc, when both the momentum and the mass term 
contribute equally to the frequency, 

fee = acMawitc) = aaH,/V2, (28) 

where Oc = a{tc) and He = H(tc). After the time tc, the sub-horizon evolution of a mode with k < kc is completely 
determined by the mass term. The second scale we define is 

ko = aoMGw{to), (29) 

that characterizes the momentum of the mode for which, the mass term starts to dominate from today on, t > tg. If 
kc < ko (or equivalently, by Eq. , tc < to), then the modes which are affected by the mass term will be within the 
observable universe today. In Fig. [TJ we show the characteristic momenta kc and fcg in relation to the evolution. 
Throughout this paper we assume that 

kc < ko, Oc < ao- (30) 

If Oc < cio is not satisfied, the effective mass would always be smaller than the Hubble expansion rate all the way 
down to the present time, and would never affect cosmological evolution of gravitational waves. We also assumed 
kc < ko for simplicity. 

With these considerations, we classify the modes depending on their momenta: 



^ We suppose that the A± domination starts at around z = 0(1) and that the evolution after that does not change gravitational wave 
spectra significantly. 
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FIG. 1: Schematic plot for the evolution of aMcw and aH. We show the critical momentum fee, defined in (|28|l . for which, 
the mass and momentum contributions are equal at the time of re-entry, while fco, defined in (|29p . is the momentum of the 
mode for which the mass term became important just today. 



• Short wavelength (fco < k): modes for which the momentum is large enough to dominate the frequency after 
horizon re-entry. We thus have tk — t^^, — and Hk — . Since the mass term never becomes 
important, the evolution is indistinguishable from its counterpart in GR. 

• Intermediate wavelength (kc < k < fco): modes for which the mass term becomes dominant after horizon re- 
entry, but before today. We still have tk — tf^, at — and Hk — H^^. We expect slight modifications with 
respect to GR signal. 

• Long wavelength (fc < kc): modes for which the mass term is dominant at horizon re-entry. Since the momentum 
is negligible throughout the evolution after horizon exit, their evolution is the same, regardless of their momenta. 
In other words, tk for these modes is roughly independent of the momentum and all modes of this type re-enter 
the horizon simultaneously. We thus have tk — tc, ak — olc and Hk — He- The largest deviation from general 
relativity is expected in this category. 

This classification is summarized in Fig. [2l 




(a)Short wavelength: fcg < fc (b)Intermcdiatc wavelength: (c)Long wavelength: k < kc 

kc < k < ko 

FIG. 2: Examples for each type of momenta based on the classification in the main text. For modes with short and intermediate 
wavelengths (fco < k and kc < k < fco), the horizon entry time is the same as in GR, with ak ~ Ok^- However, for modes with 
long wavelength (fc < kc), the momentum term never becomes dominant and the horizon entry in the massive theory occurs 
earlier than in GR, i.e. < a^^. 

In the forthcoming calculations, especially in subsection IIV El we will need an explicit expression for a^^. In order 
to determine this expression, we use the Friedmann equation for a composite fluid with nonrclativistic matter and 
radiation, that is, 
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where the subscript "eq" denotes the time of matter-radiation equahty. (See footnote [2l) Then, the horizon crossing 
time in GR can be found by solving 

^^ = :^(1 + ^)- (32) 



2(«^^)^ V «e,/ 

Depending on the dominant fluid at t = the scale factor a'jj^^ will have a different momentum dependence, 



a 



GR I %7: for fc>fc, 



k 



(33) 



where = aeqHeq is the comoving momentum of the mode which crosses the horizon at the time of equality, in 
GR. 

C. WKB solution 

As already stated in subsection IIVAI a super-horizon mode stays constant and, thus, its absolute value is given by 

|7fc|=A(fc), ioTt<tk, (34) 
where A{k) is the amplitude of the mode at the time of its generation, and for slow-roll inflation, is given by 

Here, is the value of the Hubble expansion rate at the time of horizon exit during inflation. 

Once the mode re-enters the horizon at time t = tk, it starts to oscillate with frequency u. Assuming that the 
evolution of the frequency is adiabatic, i.e. ^ , we can use the WKB approximation to obtain the solution 

7fc(i) = -^2^== exp ( -i I uadi-j] , for i > , (36) 



where C{k) is a fc-dependent integration constant. By using the thin- horizon approximation, with which we suppose 
the mode re-enters the horizon sufficiently fast, we can match the super-horizon frozen solution ([M]) and the sub- 
horizon oscillating solution l\ZQ\ at t ~ tk, resulting in 



- fort>t,, (37) 



A{k) y ij{t) a?(t) 
where loi, = i^{tk) = Hk- 

For comparison, the corresponding solution in the massless case (GR) is obtained by simply replacing — >■ fc/a and 
tk -> ifc ^ as 

D. Power spectrum 

Observationally speaking, what we are interested in is the amplitude of sub-horizon gravitational waves today as 
a function of frequency. We thus would like to compare the amplitude of sub-horizon gravitational waves in massive 
gravity with that in GR at the same frequency today. For this purpose we define the power spectrum of sub-horizon 
gravitational waves as 



r{u;o)= — (7jj7*- 

a In cjo 



'^0 -'"gw.o 



^— \jk{to)\ , k = ao^oj^ - A4iy,o : (39) 
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where 

wo = i^{to), Mawfi = MGw(to), (40) 

and we have used 

dink cjn 



din Wo ^l-^GWfi 
The power spectrum in GR with the same frequency is 



where k' corresponds to the momentum of the mode in GR with frequency loq. 
By using the WKB solution given in the previous subsection, we obtain 

k'ak\^ iOkak 



and 



ao 

where we have defined the primordial power spectrum by 



(41) 



VGR{iOo) = ^\lk' im , fc' = ao^o, (42) 



V{uj,)=['-^\ Z^Vp„mik), (43) 
' kao J ujoao 



GR N 

Pgr{ujo) = ( ^ ) Vp„m{k'), (44) 



Vp„m{k) = -^A\k). (45) 



E. Enhancement of gravitational waves and cutoff 

We are interested in the ratio of the power spectrum in massive gravity to that in GR at the same frequency today 
ujQ. From the results in the previous subsection, we obtain 



"P(^o) _ Vprtmjk) 



^Gi?X'^o) Vprira{k') 

where 



7t5^(wo), (46) 



'^(^o) ^ , GR \ ' = "oa/wo - M|^p, fc' = aoWo. (47) 

k a):, V woao V 

In the rest of this subsection, we estimate the enhancement factor 5(ajo) for the three classes of modes introduced in 
subsection IIVBI 



1. Short wavelength modes: ko <^ k 



For these small-scale gravitational waves, the frequency lj defined by (|24p is dominated by the momentum through- 
out the evolution, i.e. uj ^ k/a. Thus we have 



ujk ^ — , Wo ~ — , k' ~ k. 



Also, as illustrated in Fig. 2(a) the horizon re-entry occurs at the same time as in the GR analogue. 



„ ^ „GR ^ „GR 



Thus, this case docs not give rise to an enhancement: 



5(wo) ~ 1. 



(48) 

(49) 
(50) 
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2. Intermediate wavelength modes: kc <^ k <^ ko 
For the intermediate-scale gravitational waves, frequency at the time of horizon re-entry is dominated by momentum, 

LUk^-, (51) 

ttfc 

giving rise to the same re-entry time as in GR 

ak ^ a^. (52) 
On the other hand, the mass term comes to dominate at later time, yielding for the frequency of the modes today 

k 

u>Q ~ Mgw.o = — , k' ~ko. (53) 



Thus, these modes will undergo an enhancement: 



By using the formula ([55)) . we obtain 



'^^^^Xl^.-'^y'" k«k,,«k, . (55) 



^5/4 

1 1 for fc < fco < keq 



3. Long wavelength modes: k <^ kc 

When these modes re-enter the horizon, the mass term is already dominant in their frequency 

wfc ~ Mawitk), (56) 

and thus the re-entry time is independent of the comoving momentum. In other words, all modes with k <^ kc re-enter 
the horizon simultaneously when the Hubble expansion rate and the effective mass of gravitational waves agree 

H{tk) ~ Mawitk). (57) 

Up to factors of order one, this is essentially the definition of ic introduced in subsection IIV Bl and thus tk — tc- Hence, 
we have 

k 

ak ^ ac, u}c ^ Mcwitc) = —■ (58) 
Since we have assumed that aMcw is a growing function of time, Mqw dominates cj at t = to as well: 

Wo — Mawfl = — , fc' ~ ko. (59) 

flo 

Thus we obtain 
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4- Summary of enhancement factor 

As one can see from the schematic plot in Fig.[3l qualitative behavior of the amplification factor S{uJo) is summarized 
as follows. 

• No modification in the high frequency range (fco < k <^ ^/2AIGWfl < ^^o): S{ujo) stays almost unity. 

• Modest enhancement in the intermediate frequency range (fcc < k < ko <^ < ojo < V^Mowfl)' S{ujo) is 
proportional to some positive powers of (loq — Mg^-g)"^, and thus increases as loq decreases. 

• Sharp peak just above the cutoff (0 < fc < fcc Mgw,o < < ^c)- S{ujq) is proportional to (wg — MQ^g)"^/^, 
and thus it diverges in the limit uq — > Mgw,o- 

• No signal below the cutoff {ujq < Mcwfi)- S{ujo) = 0. 




FIG. 3: Schematic plot of the amplification factor S{u}q). 

We emphasize that the qualitative feature of the spectrum shown in Fig. [3] is quite universal; as long as aMQw is an 
increasing function of time and the two scales fcc < fco are well separated, it applies to any fiducial metric, or even 
other models of massive gravity. 

Before going into the next subsection, we briefly comment on the case where Cg is not constant."^ We assume that 
Cgk/a becomes equal to each of H and Mqw only once, similarly to the case of constant Cg we have discussed. In this 
case, Eqs. arc altered into (gHl) and into 



Cg{to)ko = aoMcwito)- 



(61) 



We redefine the horizon crossing time tk using u! in this equation. Using these new definitions, we may repeat the 
derivations in the previous subsections to obtain the formula for S{u!o). Its expression is given by 



k'ak 



Cg{to)ka'j^r V ^oao 



Cg(*fc) 



jm 



(fco < fc) 
Cs(t^)kok^ (fc < fcc 



(fcc«fc«fc-0) 



(62) 



Note that also the primordial spectrum Vprim{k) may be shifted if Cg{t) is shifted from the unity in the inflationary 
era. 



^ Note that Cg is bounded from below as Cg — 1 > —10 if we take into account the gravitational Chercnkov radiation [stI. l38l| . See 
also, e.g., Refs. [39l. |40|| for models which may give Cg 7^ 1. 
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F. Signatures in the spectrum 

111 the previous subsection, we have shown that the enhancement factor S{ujq) defined by (|46p diverges in the hmit 
luq — > Mqiyq + 0. However, the power spectrum V{ll!q) itself docs not diverge, provided that 

lim k^'^Vnrimik) < oo. (63) 
fc->+0 \ / 

This is the case, for example, if the primordial gravitational waves are generated by inflation in the early universe 
and if the e- folding number of the inflation is finite. Nonetheless, from the divergence of 5(cjo) it is likely that the 
power spectrum 7'(wo) has a sharp peak just above the cutoff at ujq = Mcwfi- Below the cutoff, the power spectrum 
should vanish. 

The peak just above ujq — Mgw.o may become rather high if the primordial spectrum Vprim{k) of gravitational 
waves continues to a far infrared. For example, this situation can be realized when the duration of inflation is 
prolonged. Concretely, Vprimik) is expected to continue down to fc ~ aoi/oe"^""'''''' , where N^xtra is the number of 
additional e-folds to the minimum needed for the spatial curvature to stay below the observational bounds. Hence, 
by the formula (j60p . the enhancement factor at the peak is estimated as 



Future observation of gravitational waves may reveal deviation from the spectrum predicted by GR. If the cutoff 
and the sharp peak just above it are found, then this might be a signature of massive gravity. In this case, from the 
frequency at the cutoff one can read off the value of Adcwfl- It might also be possible to read off some information 
about the fall-off behavior of the primordial spectrum in the far infrared (j63p . from the height of the peak. 



V. NUMERICAL RESULTS 



In the previous section we have shown that the power spectrum of stochastic gravitational waves in massive gravity 
vanishes at frequencies below a cutoff and a sharp peak is expected just above the cutoff. The cutoff frequency agrees 
with the effective mass of gravitational waves today Mawfi- It is therefore interesting if the effective mass M^vi/.o is 
within the frequency range to which gravitational wave detectors are sensitive. For example, space-based gravitational 
observatories like the Laser Interferometer Space Antenna (LISA) [4l[ is sensitive to the frequency range 10~^ — 1 Hz. 
The DECi-hertz Interferometer Gravitational wave Observatory (DECIGO) [i^l has better sensitivity in the frequency 
range 10~^ — 1 Hz. For lower frequencies, pulsar timing arrays such as the Square Kilometre Array (SKA) [43] and 
the Parkes Pulsar Timing Array (PPTA) [4J| have sensitivities down to 10^^ Hz. The combined sensitivity range for 
these observatories is thus 10^^ — 1 Hz. Thus, using these observatories, one might be able to hope future detection 
of the cutoff and the peak in the gravitational wave spectrum if 

10"** Hz < Mawfi < 1 Hz. (65) 

The spectrum obtained analytically in the previous section, especially its peak structure, is not sensitive to the history 
of Mcwit) in the early universe but is essentially determined by its present value Mgw.o- Thus, for the purpose of 
numerical calculations, we consider time-independent Mqw = ^GW,o in the range ([GSp ."^ 

In Fig. |4l we compare the numerically obtained gravitational wave spectrum and the analytic estimation derived 
in the previous section. The former is obtained by numerically solving Eq. ^ with the mass term assumed to be 
Mcwfl = lO^^iTo ^ 2 X 10^'^ Hz, for which we find fco/^^o = 10^° and kc/Ho ^ lO'*. The evolution of the gravitational 
wave inside the Hubble or mass horizon is calculated using the WKB solution, Eq. (|36p , in order to save computational 
time. The latter is obtained by multiplying the numerically obtained spectrum in GR by the analytically obtained 
factor 5(ajo), given in Eq. In both cases, we assume that the evolution of the Hubble expansion rate is determined 
by the density parameters of radiation Vlrh? = 4.15 x 10~^, matter ilmh^ = 0.1344, and cosmological constant (either 
the genuine cosmological constant or the effective one induced by the graviton mass term, or their sum) f^A = 0.734, 



Note that in Sec. |ll] we have assumed that propagating degrees of freedom in the gravity sector are the two tensor modes only. This 
means that the scalar and vector parts of linear perturbations behave exactly like GR and thus there is no constraint from modification 
of Newtonian potential as far as linear perturbation is a good approximation. Sec e.g. Rcfs. [45l for the discussions and bound on 
the graviton mass when those degrees of freedom arc taken into account. 
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with the normalized Hubble rate h = 0.71 [47|. The primordial spectrum is assumed to be scale-invariant and its 
amplitude is normalized by 'Pprim(fc) = 2.43 x 10~^°, where we assumed the tensor-to-scalar ratio to be 0.1 and the 
amplitude of the scalar primordial spectrum to be 2.43 x 10^^ [13]. We also suppose that the primordial spectrum 
have an IR cutoff at k/{aoHa) = e"^=-"-" = 10"^ We can clearly see that the numerical result matches almost exactly 
the analytic result derived in Sec. IIVI and the peak height is well described by Eq. ([M]) . which gives Speak ~ 10^^. 



10-^ 



10-' 



numerical 
semi- analytic 



10"^ 



10"' 



10" 



10'' 



10" 



10" 



10" 



10-' 



10" 



COq [Hz] 



FIG. 4: Gravitational power spectrum in the massive gravity with Maw,o = W^" Ho ~ 2 x 10~^ Hz shown with respect to ujo- 
The solid (red) curve is the spectrum obtained by the numerical calculation, and the dashed (blue) curve is the semi-analytical 
spectrum given by Eq. (|46|l . where the GR spectrum is calculated numerically and the enhancement factor 5((xjo) is given 
analytically by Eqs. ^0^, and ((60)) . 

Since the plot in Fig. |l]is shown with respect to wq ~ [(^/"o)^ + ^gwoV^"^' most of the peak structure in the 
spectrum is compressed to a narrow region at ujq ~ Mgw,o and its detailed structure is too subtle to see. To make the 
fine structure of the peak manifest, we show the plots with respect to k in Fig. [5] We can see that the k dependence 
changes at fco and kc as prescribed in Sec. IIVI and that the semi-analytical result well reproduces the numerical one. 




k/a(,Ho 



FIG. 5: The gravitational power spectra in the massive gravity theory with Mgw,o ~ lO^^-ffo ~ 2 x 10^^ Hz shown with 
respect to k/aoHo- The spectra are obtained by the numerical calculation (solid red) and semi-analytic estimate by Eqs. (|46p . 
(|50|l . (|55|) and (|60|l (dotted blue), respectively. For the latter, the GR spectrum is calculated numerically and the enhancement 
factor S{uJo) is calculated analytically. 

Before closing this section, we briefly comment on the observability of the peak feature in the spectrum. As seen 
from Fig. [HI the peak is observable by the pulsar timing or eLISA, and/or DECIGO if the mass is in the range of 
Eq. (p5)) . For example, for the case of Mgw.o = lO^^Hz and a primordial spectrum with an IR cutoff at fc/oo = ^/Hq, 
we find that the enhancement factor S{u!o) reaches ^ 10^^. Since the enhancement is so large, we can expect the 
peak to be observable even when the tensor-to-scalar ratio takes much smaller value. If Mgw,o is much larger than 
1 Hz, the peak in the spectrum would be more manifest, and the signal would be observable even for ground-based 
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FIG. 6: Gravitational power spectra for Mcwfi/^TT = 10^®, 10^^, 10~*Hz shown with sensitivity curves of future experiments. 
The curves are plotted assuming 10-year runs of pulsar timing array with SKA [i^, the single detector sensitivity of eLISA 
[49I I. and the cross-correlation analysis of a 3-year observation with DECIGO For the primordial spectra, we consider the 
tensor-to-scalar ratio to be 0.1 and e~^=°=*'"'' = 10"^. 

observatories, such as LIGO or KAGRA [50l.[5l|. 

Note that, due to the finiteness of the observation time Tobs, the observed peak in the Fourier space is smoothed 
with a width of A/ ~ 1/Tobs- The signal would be detectable even after taking this effect into account if the 
amplification is sufficiently large. For example, following the calculations in Ref. [s^ which discussed the pulsar timing 
observations, we find that the total RMS fluctuation of the timing residuals induced by the stochastic gravitational 
wave background is amplified relative to GR by a factor ^ 10"^ x {Tobsl Ho)~^{Maw,ol Ho)~^ \og{e^'"^*'''^ Mow,o I H^) ■ 
For Mcwfi ^ 10~^ Hz, this factor becomes large enough for the signal to be detectable by the cross correlation 
analysis using PPTA argued in Ref. |52| unless the tensor-to-scalar ratio is extremely small. 

VI. SUMMARY AND DISCUSSIONS 

In this paper, we discussed how to probe the graviton mass using the gravitational wave observations. Our results 
are based on a general quadratic action given by Eq. ([S]) in which the graviton mass can be time-dependent, and thus 
they apply universally to any massive gravity theories in which modification from GR appears primarily in the tensor 
modes. The salient feature of the result is the sharp peak in the gravitational wave spectrum, whose position and 
height may tell us about the graviton mass of today and the length of the inflationary period. In Sec. |Vl we conflrmed 
our analytical results match well with the numerical results which do not assume the thin-horizon approximation, 
and also argued the observability of the gravitational wave signal for various observatories. It would be interesting 
to apply the results to explicit examples of nonlinear completion of massive gravity theories and to derive constraints 
on such models based on the gravitational wave observations. 

Although the example discussed in Sec. IIIII (as well as in Ref. [s^) was recently shown to suffer from nonlinear 
instabilities [3l| . it is expected that other working examples exist. A construction with no extra polarization with 
respect to GR was discussed in (soj . where a constant mass contribution appears in the dispersion relations of 
the gravitational waves. Another potential example is the scenario considered in (35| . where the scalar mode is 
rendered redundant by introduction of some symmetry. Although the vector graviton degrees are still present in this 
construction, the lack of the scalar graviton would relax the bounds from solar system tests. 

Below, we briefly discuss other probes which are sensitive to MQw{t) at the time other than today. One interesting 
probe is the cosmic microwave background (CMB), especially its B-modc polarization spectrum which is sourccd 
solely by gravitational waves. The signatures in CMB spectrum in the case of the constant mass Mqw was discussed 
by Ref. [53| • Following the argument to derive the B-mode spectrum and generalizing it to the case of time-dependent 
mass terixL, we find that the main contribution to the spectrum is sensitive to the value of the mass at the time of the 
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recombination^. The most interesting observational feature will be obtained when the graviton mass at that time is 
in the range (lOMpc)"^ < Mqw < (lOKpc)"^. For such a value of Mqw: the contribution to the B-mode spectrum 
from the gravitational wave increases compared to the GR case for the modes of low angular multipole i, and a plateau 
in the spectrum up to £ < 10""^ x Mow/Hq ~ 100 is expected to emerge. There will be no signal for smaller values of 
Mqw, while for larger values, the entire spectrum will be suppressed, since the rapid oscillation of the gravitational 
wave, driven by the large Mow, will average out the contribution to the CMB polarization spectrum. 

Another interesting probe of the modification in tensor sector of gravity is the inspiralling compact binaries; the 
changes in the propagation speed due to the mass term will give rise to a modification in the evolution of the phase 
of the gravitational waves emitted by those binaries. Such an effect and the constraint on the constant graviton mass 
are argued in Ref. [s^l • It may be interesting to generalize this analysis to the case of time-dependent graviton mass 
and discuss its observability for the future-planned gravitational wave observatories. 

Finally, we give a brief comment on the time evolution of the gravitational wave energy density and a possible 
constraint on it. The energy spectrum of the gravitational wave defined with respect to wavenumber k is given by 



where w is given by Eq. p4| and the final expression is valid for the WKB solution of Eq. ((36)) . The ratio of the 
energy spectrum of the massive gravity relative to GR can be obtained following the derivations similar to Sec. lIVI To 
discuss the constraint on the gravitational wave energy density coming from the Big Bang nucleosynthesis (BBN) , we 
need to do the analysis taking t to be the time of BBN. We assume that the gravitational wave primordial amplitude 
is given by Eq. ([55)1 . and consider the case that MQiY{t) becomes larger than the Hubble scale before BBN. In such 
a case, the ratio is given by 



where fco and kc are defined in a similar manner with Sec. llVl replacing to by teBN- Roughly speaking, the total energy 
density of the gravitational wave becomes larger than that in GR by a factor of k^/kc = Mq\y (<bbn ) / \/ Hp,PMMqw{tr) 
if kc is sufficiently large. Since the total energy density must satisfy J d{\ogk)^Gw{k) < 10^^ at least [55[, this 
enhancement of the gravitational wave energy density may give a tighter constraint on the primordial amplitude of 
the gravitational wave or the inflationary models from which it is derived. Obviously this constraint will be sensitive 
to McwitBBN)- It may be interesting to pursue this issue using explicit examples of the massive gravity theories. 
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In the nonlinear theory of massive gravity, the mass term not only introduces the effective mass to gravitational 
waves but also drives self-accelerating universes. In order to attribute the late-time acceleration of our universe to 



^ We find the contribution of the gravitational wave to the CMB B-mode spectrum is roughly proportional to (fc^ + Mqw (v)) drj, 
where rjr is the conformal time at the time of the recombination. In a generic case, this value is sensitive to MowiVr) rather than 
^Gwiv) ^t any other r}. 




(66) 




(67) 
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Appendix A: Example from nonlinear massive gravity 
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the self- acceleration from massive gravity, the effective cosmological constant A± due to the mass term should be set 
to the observed value: 

A± ^ Hi (Al) 

Therefore, a relevant question now is whether the nonlinear theory of massive gravity can simultaneously accommodate 
the effective mass in the preferred range ((65)) and the preferred value of the effective cosmological constant (|Aip . These 
two quantities are given by the formulas ([^ and (PT|) , respectively. (For the sign of the effective cosmological constant 
in each branch, see Figure 1 of Ref. jssj.) Thus, for generic values of the parameters and Q!4, the graviton mass 
nig in the action should be 

nig - Ho, (A2) 
and the quantity r = H/{X±Hf) at present must be in the range 

10^° <ro< 10^^ (A3) 

where ro = r{to). 

If we suppose that the primordial spectrum Vprim{k) of gravitational waves is generated by inflation then in order 
for the inflationary mechanism of generation of quantum fluctuations to work, the effective mass of gravitational waves 
must be sufficiently lower than the Hubble expansion rate during inflation. This is translated to the condition 

H^ 

r <^ — |- during inflation, (A4) 

^0 

where Hj is the Hubble expansion rate during inflation. 

For example, if we demand the Poincare symmetry in the field space as a global symmetry of the theory then 
the fiducial metric is restricted to Minkowski and the FRW background of the physical metric is inevitably an open 
universe [3l ■ In this case we have 

(A5) 



where K (< 0) is the spatial curvature constant. It is evident that r increases during inflation and decreases after that 
(see footnote [2]). Hence, the condition (|A3P implies that r ^ 1 throughout the evolution after inflation. In this case 
^GW ~ rrri^ for the '+' branch (or Mqt^ ^ —rrrig for the '-' branch), and the assumption is actually satisfled, 
provided that P < p and > (or < 0, respectively). 

Since r decreases after inflation, the condition (|A4[) leads to an upper bound on rg. We thus need to see if such 
an upper bound on tq is consistent with the preferred range (jASp . For this purpose, supposing that the universe is 
dominated by inflaton oscillation (with the effective equation of state parameter w = Q) between the end of inflation 
(a = a/) and the onset of the radiation dominated epoch (a = a^), the scale factor today is estimated by the entropy 
conservation as 

^ fanV'' ^ (A6) 

a/ To V «■/ / 

where Tg is the photon temperature today and g is the number of relativistic degrees of freedom that eventually inject 
entropy to photons. By using this formula, the bound (|A4|) at the end of inflation t = tf is rewritten as 



Hq ao Hj ao 1/12 V-^-^P'^-f f o-R 



1/4 



^o=rf— <— — — , (A7) 

Hi af Ho af To Ho V «/ / 

where Vf = r(tf). For example, for g^^^'^ = 0{1) with the instantaneous reheating approximation a^ ~ a/, this leads 
to 



and is consistent with the preferred range (jA3[) . Thus, in this simple example with Minkowski flducial metric, the 
range (jA3[) is consistent with the inflationary generation mechanism of primordial gravitational waves. 
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Note that, by the formula (jASp . the extreme hugeness of tq imphcd by the preferred range (|A3I) is translated to 
the extreme flatness of our universe today: 10^^^ < \ fttot ~ 1| < 10~^°. Thus, the number of e-folds during inflation 
must be much more than what is needed to solve the flatness problem. The preferred range (jA3[) corresponds to the 
number of extra e-folds ~ In in the range 

45 < Ne^tra < 65. (A9) 

With this large number of extra e-folds, the primordial spectrum Vprimik) of gravitational waves is expected to 
continue to a far infrared. Therefore, from the discussion in subsection IIVFI we conclude that the peak just above 
Wo = Mcw,o should be rather high. Concretely, Vprim{k) is expected to continue down to fc ~ aoi?oe~^°^''^° . By the 
formula (|60l). this corresponds to the large enhancement at the peak. 



'-'peak ~ r'R TT^^ rn TT^'O- {■^^'^J 



Appendix B: General fiducial metric 

1. FRW-type fiducial metric 

The background equation of motion for a general fiducial metric is derived in the Appendix A.l of Ref. [11]. We 
briefly summarize that derivation for reference. 
We consider a FRW-typc fiducial metric given by 

- ^h\^°)d^^°dj° + d2(0O)f)(f )a^0»a.<^^ (Bi) 

where n and a are general functions of 0°. Applying the field redefinition given by 0') = (/((/)*'),(/)'), where / is 
a function to be determined, and taking the gauge 0^ = .t^, the fiducial metric of Eq. (|Bip becomes 

U,dx^'dx'' = -h^{f{t))f{tfdt^ + a^{f(t))n^pdx'dx^, (B2) 

where the dot (*) denotes the derivative with respect to t. For this fidueial metric, the matrix /C^^ becomes 

/C^„ = Diag(C,^,^,0, C = l--^, e = (B3) 



Ig = / d^xVU (sKNa - + mlLg ] , (B4) 



Then, up to boundary terms the gravity action reduces to 
where 

Lg=Na^{C2+a3C3 + a4C4), C2=m + 0, C3 = e{K + 0, A = (B5) 
The variation of this action with respect to / yields 

= + _ f ^ 3Na^ (aH aHf) [l + 2^ + a,{2 + + «4C'] , (B6) 



ac df df V dQ df 

where H and Hf are the Hubble parameters of the physical and fiducial metrics defined by 

H{t) = 4r^ Hf{f{t)) = ^ (a' = ^\. (B7) 

Below, we focus on the solution of Eq. ()B6P given by^ 

l + 2e + a3e(2 + + a4e' = 0, (B8) 



See Refs. 1341. ISal for cosmological solutions on the other branch. 



18 



whose solution is given by 



KM 

ait) ' 



_ 1 + 2a3 + Q!4 ± ^1 + as 



a4 



as + ^4 



(B9) 



Varying the action with respect to N and a and using the solution of Eq. (|B9p , we obtain Eq. ^U\\ as the background 
equations. The second equation in Eq. (|20p is consistent with the first one if the matter fluid obeys usual conservation 
law p + 3{a/a){p + p) = 0. The evolution of a is determined by Eq. ((20)l for given N, and it determines f{t) via 
Eq. (jB9[) . Note that n(/(t)), a{f{t)) and Hf{f(t)) correspond to n(t), and of SecHIl respectively. (Almost 

equivalently, we may set /(t) = t and regard Q;(i) and n(i) are given functions. In this case, Eq. ()B9[) and the 
Friedmann equation in Eq. arc regarded as equations to fix a{t) and N{t), respectively.) 

One example of the fiducial metric of this type is the Minkowski metric in the open chart (K < 0) given by [isj 



where we used the field redefinition given by 



(BIO) 



(Bll) 



Taking the gauge <jf = , we find Eq. (jB9[) becomes X± = ^J\K\f{t)/a{t), which determines f{t) in terms of a{t). 
The Friedmann equation, Eq. ((20|) . holds as it is. In this case we have n(0°) = \f'{4>^)\ and a{(lP) ~ yJ\K\f{4f'), and 
thus Hf = sign(/')// = sign(/')v^/^±a- 



2. (Anti-)de Sitter fiducial metric 

In this section, we comment that the fiducial metric of Eq. ([T7| for any K can be obtained from a single dc Sitter 
metric, and the for K > can be obtained from the anti-de Sitter metric. 
We begin with the de Sitter fiducial metric with fiat spatial section given by 



(B12) 



By the redefinition (p'^ = /(</)") and = this fiducial metric takes the form of Eq. ([T7)) with K = 0, Hf — Hf, 
n = |/| and a ^ exp(2iJ//(0O)) . 

Assuming K > 0, the field redefinition given by 



sinh/(<^") + coshf{(j)")^l~ K5^j(j)^(l)3 



tanh/(^o) + ^1 - KS,j(f>^(j)^ 



makes the fiducial metric of Eq. (|B12p into 



Uu = HJ^ (-f'\^°)d^^°d,cb° + A-coshV(0°) n[pd^^'d,(ly' 



(B13) 



(B14) 



which is of the form of Eq. ^7} with AT > and Hf = Hf tanh/((/)°). The counterpart of Eq. ((B9| for Eq. (|B14p is 
given by 



and the Friedmann equation, Eq. ((20)) . holds as it is. 
Similarly, for K < 0, the redefinition given by 



i^cosh(/(i)) 



Hf 



ait) 



(B15) 



makes Eq. (|B12I) into 



cosh/(0°) + sinh/(/)y^l - K5,j4)'(t>> 



\Kl 



coth/(^o) + ^/i^na~^ 

.U. = HJ^ [-/'' (00)9^009,00 + |X| sinh2/(0°) f^lf ^9^0^9,0^' 



(B16) 



(B17) 
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which is of the form of Eq. (IT7|) with K < and Hf = Hf coth /((/>'^). The counterpart of Eq. (jB9[) is given by 

x,.mi^. (BIS, 

Hf a[t) 

Next, we consider the anti-de Sitter (AdS) fiducial metric given by 

= - (1 + fc'^'^V-'') <9^¥'"<9,</." + ^7|7'■')9^(^■'9,v'^ (B19) 
which is the AdS metric with curvature k in the global chart. Assuming K < 0, the field redefinition given by 

ip° = k arcsin ( '^"^^'^"^ | , = cos/ (0° ) 0' (B20) 

make the fiducial metric into 

- + |X| cos2/(</'") f^lf , (B21) 

for which Hf = fctan/((/)°) and the counterpart of Eq. (|B9|) is 

k a[t) 

We summarize these results in Table HI Note the the constraint equations, which are the counterparts of Eq. (jB9p . 
give restrictions on the trajectory of a{t). For example, the constraint equation in the K > case allows only the 
trajectory of a{t) which starts from the infinity, decreases to the minimum and grows to the infinity again. Even 
though the a(t) which starts from zero, increases to the maximum and comes back to zero is also allowed by the 
Fricdmann equation, it is not consistent with the constraint equation and is not a cosmological solution in this model. 



TABLE I: Results summary. 











Hf 


Constraint eq. 


dS 


K>Q 


Hi' log 


sinh/(/) + cosh f{(l)°)-s/l - KSij(l)icj>^j 


tanh/(0O) + ^l-A'4,j0i<#.J 


/f/tanh(/(/)) 


Y -sfK cosh /(t) 
Hj a{t) 


K = Q 




0' 


Hf 


^ oxp(%/(t)) 
^± - a(t) 


K <0 


Hj'log 


cosh /(/) + sinh f((l)°)^l - KSijfjy^ij)^ 


coth f{<l>'>) + y/l-KSij4,^<l>i 




Y ^/\K\ sinhf(t) 
Hf a(t) 


AdS 


K < 


k arcsin 


( Auf{4>°) \ 




^cos/(0«)<^» 


fctan/(<j!>°) 


Y, _ cos/{t) 






^± k a{t) 



Before closing this appendix, we briefly mention some examples of Mqw- For instance, in the case of the + blanch 
solution with the Minkowski fiducial metric, 

^GW I Minkowski °^ (° + constant) . (B23) 

As a result, a^|Mg.^y| always grows with time, provided that the dominant fiuid's equation of state is w < 1. Another 
example is the de Sitter fiducial metric, for which we have 



^'GwLr.^;ur.r « ( - + constant ) . (B24) 



a 



In this case, the contribution from the first term in o^Mq,^ will grow if w < 1/3. Although for radiation dominated 
universe, this contribution is constant, the second term will still grow as a^. In the light of these examples and the 
late-time stability of the system, we adopt the assumption ([25]) throughout this paper. 
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Appendix C: Analytic solutions of 7^ 

We show the analytic solution for Eq. ([9]) for power-law type Mqt^ = A/^q(a/acq)^, which covers the simple examples 
such as the constant and also the examples mentioned above. The aim of this section is to confirm the validity 

of the thin-horizon approximation used in Sec. IIVI 

We define the conformal time -q from the scale factor a ~ acq {t/Uc^f^ (RD) and a — o,oq{t/tcq)^^^ (MD) as 

7o ait) ||(_^) _i (MD) acq ^ ^ " ■ — 

where a is a constant and Tyeq = 2ateq/acq- 

We consider the case that Mqw becomes dominant in the MD era first. Eq. ^ in this case becomes 




7fc 



where we introduced k — kij^q, M^q = aoq'7cq-^cq, o, = a/ acq and fj = — V \ = 2a^/^/3. Below, we use an 
approximation depending on which of k and aMcw is dominant at the moment of the horizon crossing. 

1. fc 3> aMcw at the horizon crossing: 

In this case, Eq. (|C2p is approximated near the moment of the horizon crossing as 

^H-(^^-^)% = 0, (C3) 

which is same as that in the pure GR. This equation has an exact solution which becomes regular for ki) 
given by 

7. ^ C X I cos(i.^) + ^ Cf. (C4) 

The primordial amplitude is given by \^k\ ~ k~^/^(H^/AIpi) and 

I- I H^a 4iJ*acq77cq^ .2 _ - -2 /r^c-N 

and then wc may identify it with the right-hand side of Eq. (|C4p in the limit of kfj <C 1. The solution in the 
WKB regime, on the other hand, is given by 

Ik = cos ( / kdrj] , (C6) 



where k^ = (k^ + Mq^^)'^^^. Comparing it with Eq. (jC4p in the limit of ki) ^ 1, we may fix the coefficient C" 
to have 

7fc = ^77: cos ( / kdfj] . (C7) 



fc3/2 

Let us compare Eq. (|C7P with 7^ obtained using the thin-horizon approximation. The condition to determine 
a at the horizon crossing is given by fc = aH = a' /a = 2/77. Matching the WKB solution with the primordial 
solution at rj defined by this equation, wc find the solution is given by 

which coincide with Eq. (|C7p up to a numerical factor. 
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2. aMow ^ k at the horizon crossing: 
Eq. (|C2[) in this case is approximated as 



which have an exact solution that becomes regular for fc?} — > given by 

(|)^-tglr(i + |) /3(|)^„A ,^0 

7fc = C X ^ t' J3 -rj^- C-q-^, (CIO) 

A/el n « / 

where J is the Bessel function and we have defined q = 2{p + 3). Fixing C using Eq. (jCSp . we find 7fc to be 

2 g „77F / 1 _L 3 



7fc = 7* X ^ rf/Va ^j}^ 

Afel ^\ 1 J 

,^oo g'^^r(l+|) ^.(D^-tjOeqi /3(|)^.A ^^^^^ 
^= ^ X COS — l^ii ^7^2 ^ (Cll) 

where we neglected the phase of the oscillation. This expression is valid as long as Mow dominates over k/a. 

Below, we compare it with the expression obtained from the thin-horizon approximation. The condition to 
determine rj at the horizon crossing is given by 

'2\''^+'^ ~ . ^ 2 1 /3\Si „ 1 

I) M^^f(^+'^ = ^ .-. ,, = 2^(^0 Af-q^. (C12) 

Matching the WKB solution with the primordial solution at this 7], we find 

3_3 , _3 

_^ ^ X 7. (i)- -Afcq- / /• I \ ^^^^^ 



This expression coincides with Eq. (jCll| up to an 0(1) factor for generic q > 0, and thus the thin- horizon 
approximation is fairly good even in this case. 

Finally, we briefiy comment on the case that Mew becomes dominant in the RD era. In this case, Eq. © becomes 
7" -t- (fc^ -t- a^A/g^) 7 = 0. As long as fc^ -I- a'^M^^^ is positive and its time variation is sufficiently slow, the solution 
will be well approximated by the WKB solution. The amplitude of this solution is determined by matching it to the 
primordial solution. By construction, this solution will be well approximated by the solution given by the thin-horizon 
approximation in general. 
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